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Abstract 

The Lenz-Ising-Onsager (LIO) problem in an external magnetic field in the 
second quantization representation is the subject of consideration of the paper. 
It is shown that the operator Vh in the second quantization representation 
corresponding to Ising spins interaction with the external magnetic field H 
can be represented in terms of single-subscript creation and anihilation Fermi 
operators in such a form that the operator Vh commutes with the operator P = 
(—1)'^, where S = J2mPmPm is the operator of a total number of Fermions. 
The possible consequences of such representation with it's relation to the LIO 
is discussed. In particular, the constructive proof of the Lee- Yang theorem 
on the absence of phase transition for Ising model in nonzero magnetic field 
7^ 0) is demonstrated. 

PACS number(s): 05.50.+q 
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I. INTRODUCTION 



As is known, up to now there is no solution for the 2D Ising model (IM) in an external 
magnetic field H (Lenz-Ising-Onsager problem - LIO) still. The algebraic, topological and 
graph-theoretical difficulties of taking into account the external magnetic field in LIO prob- 
lem are well known and widely discussed |]T]-0- As is known |]^, the statistical sum for 
the 2D Ising model in external field {H) in the representation of second quantization can 
be written in the form: 

Z2D{h) = Tr{Vr = Tr{V,V,V^r, (1) 

where the operators Vj, expressed in terms of the Fermi creation and annihilation operators 
(c]„. Cm), are of the form: 

= (2sinh2iri)*^/2exp [-2K* =l(4^c„^ - 1/2)], 
V2 = exp {K2 [E^lti (4 - Cm){cUi + Wi) - (-1)*^(4 - cm){c\ + ci)] } , (2) 
Vh = exp {hJ2m=i exp [iir E^=i 4cp\ (cj„ + c„) } , 

where Kj = jSJj] (j = 1,2,) and h = jSH (/? = l/ksT - inverse temperature, H - vnexnee 
magnitnoe pole), and M = cl^Cm is the operator of the total number of particles and 
and Ki are connected by the following formulae: 

tanh(/s:i) = exp(-2is:*), or sinh 27^1 sinh 2/s:* = 1. (3) 

One can see that the operator Vh in the second quantization representation, that describes 
interaction of the spins with external magnetic field, has rather complicated structure. It is 
easy to see that this operator does not commute with the operator P = (—1)^. As a result 
the operator V2 has also not a very tractable form, i.e. it has not the needed translational 
symmetry (2). More exactly, although the operators Vi and V2 commute with the operator 
P, the operator V (1) does not commute with the operator P, i.e. [P, V]_ 7^ 0, because 
[P, Vft,]_ 7^ 0. Therefore, we can not divide all states of the operator V = ViV2Vh into 
eigenstates of the operator P with eigenvalues A = ±1, and this leads to nonconservation 
of the states with even and odd numbers of fermions (for details see [^). Namely this is 
the fundamental reason which stops solving the problem under consideration within this 
formalism. 

Coming back to the difficulties mentioned above which are connected with the operator 
Vh, (2), it is now clear that to overcome the troubles within the approach 0, one should 
find an appropriate method of substituting the operator Vh (2) with another one which 
would be equivalent to the former in the sense of correct counting of the interaction of 
external magnetic field with the spins of the system. Namely, as it could be easily seen, the 
only contribution to Z2D (1) from the operator Vh comes, in the representation of second 
quantization, from the "even" part with respect to operators cl^,Cm of the operator Vh- 

It was shown in the author's paper ||^ how this difficulty can be overcome by transition 
to the space of higher dimensionality and supposing then one of the interaction constants 
to be equal to zero. Then all the operators Vj in the second quantization representation are 
expressed in terms of two-subscript creation and anihilation Fermi operators in the finite- 
dimensional Fock space 2^^"^^ in the particle number representation. These operators Vj 
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commute with the operator P = (—1)"^, where S = Y^nm^Aim^nm is the operator of total 
number of fermions. Here we discuss briefly how the analogous result could be achieved in 
the space of single-subscript Fermi operators, that is, in the finite-dimensional Fock space 
2^. In particular, I'm intend to sketch briefly that the partition function for the 2D Ising 
model in an external magnetic field can be expressed in the form: 



Z2D{h) = Tr{V) 
(2sinh2Ki)^^/2 



N 



Tr{ViV2VH) 
V, = exp {K2 [E^'=i(/5l - Pm){PUl + Pm+l] 



exp 



1/2) 



(4) 



Vu = (cosh h)'' n n 1 + tanh^ h{f3l - Pm){PUk + ^^^^ 

m=l k=l 



(5) 



It is quite obvious that the operator Vh (5) commutes with the operator P = (—1)"^, where 
S = Y.m=i PpPp operator of total number of fermions, and hence, the operator V2 (4) 

also can be expressed in that form, just as it is done here, that is without term ~ (— 1)'^ 
unlike the representation (2). 



II. THE PARTITION FUNCTION 

The main idea of representation (4)-(5) is to consider the 2D Ising model in an external 
magnetic field in terms of second quantization representation, in in the space of two-subscript 
creation and anihilation Fermi operators and after that to factorize the corresponding oper- 
ators T12 (see below). Namely, it was shown in the author's that the partition function 
of 2D IM in an external magnetic field is of the form: 



zMh) = {2 cosh hf \ n Aii<i)][ n Ai{p,h)]{omhmo), 



(6) 



where 



exp 
exp 



L^n=l ^1=1 ^m=l "'K'' J HnmHn+l,m 



■'171= 

M-m ■ 



a(0 = i Eo«?<^ 25i(g) sin(/g); 
, b{k) = -ti Eo<p<n 2S2(p) sin(A;p), 



(7) 



and 



B2{P) 



A2{p, h) = cosh 2K2 — sinh 2K2 cosp + a{h) sinh 2X2 sinp, 
v4i(g) = cosh 27^1 - sinh 2/^1 cosg, a{h) = tanh^ h^^^, 

a(/i)[cosh 2_ft'2+sinh 2_ftr2 cosp]+sinh2A'2 sinp r) r \ sinh2A'isinq 

A2ip,h) ' -°HyJ - A^) ■ 



(8) 



Here (ajj^,a„m) and {Pimy Pnm) are the two-subscript Fermi operators, |0) is the fermionic 
vacuum function in the finite-dimensional Fock space of 2^^^ dimensions in the occupation 
number representation. Fermi operators a- and (3- are related to each other by canonical 



unitary transformations {a. 



nmP-nm 
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Z^fc=n+1 ^p=l 2^k=l l^p=m+l\ ^kp^kp — [• • ■\Hkph'kp- 

In order to transform the representation (6)-(8) into (4)-(5), one should make the following: 
to rewrite the operator Ti (7) in terms of Pauli operators (r^), to apply the well-known 
representation 



1 

exp(C^) = — / expK' + 2C^]de, 



(which is correct for the totally bounded operators), to turn Pauli operators (r^) into to 
the Fermi ones (aj^^,a„m), then to "drag" the emerging phase factors through |0), 0, to 
introduce Ising-type variables (/inm = il); to apply for the second time transfer-matrix 
method and finally to use the one-dimensional Jordan- Wigner transformation. Omitting 
the great number of of intermediate calculations, for the partition function (6) we have: 



ZMh) = (cosh/i)^^^ ■ Tr{VrV2V:)^, (10) 

where operators Vi^2 is defined by formulae (4) and trace is calculated by means of (3- 
operators, where operator is of the form: 

M M- ^ 

V* = {Ole'^-^'^E^^.E^Sr-^^^-- Y[[i + (al - a^)(aL+i + «™+i) ■ iPl - Pm)iPL+i + Pm+i)] |0)„ 



m=l 



(11) 



Here the single-subscript a- and /?- Fermi operators are absolutely independent ones (as they 
are introduced), commuting or anti-commuting with each other (the final result does not 
depend on the fact, as it should be expected) and the vacuum martix element is calculated 
by means of a- operators. Using Wick theorem, one can easily calculate the vacuum matrix 
element in (11) and obtain the final expression (5) for the operator Vh, taking into account 
the factor (cosh h)^^^ at the Tr{...) in (10). One can say, using the figure of speech that the 
operator exp[tanh^ hY,m=i J2k^Ji^ c^m+kOim] within the brackets (0|...|0) acts on the operator 



rimf---] just like the "Ockham razor" does, cutting off, due to Wick theorem |]T0[, all the 
superfluous terms. 

From (4)- (5) one can derive some consequences, in particular, the constructive proof 
of well-known Lee- Yang theorem |Tl|,[12[ on the absence of phase transition for the Ising 
model in nonzero magnetic field ^ 0. Indeed, it can be easily shown from (4)-(5) that 
at /i = we have Onsager solution and supposing Ki to be zero {Ki = 0), we have next 
representation for the partition function: 

Z,o = {2coshh)''{0\{V2V^)\0), 
= exp [tanh^ /^E^ii Hfj^ PI(^1+p_ 

from which we can derive exactly the classic Ising result 0. 

But it is not trivial to have by means of (4)- (5) Ising solution supposing K2 to be zero 
{K2 = 0). For this case operator V2{K2 = 0) = I and the problem is reduced to the 
calculation of trace with the operator (ViV^)^, which is not easy. However, as it follows 
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from (4)-(5) by it's construction, we have here the Ising solution also. The last one allows 



to prove constructively the Lee- Yang theorem for Ising model [TT|,[T^ . 

First of all, it is easy to demonstrate that at small but finite external magnetic field 
(/i ~ £ <^ 1), the expression for operator V/j, (5) can be represented as: 



Vh = exp 



M M-m 



m=l k=l 
-4 



(12) 



with the accuracy up to ~ e . It is obvious that the next Hamiltonian: 



N,M 



n 



n=l,m=l 



N,M M-m 



h T 

'^B-i- n=l,m=l k=l 



corresponds to the partition function (4) with the operator Vh (12), with the accuracy up 
to unessential constant (~ H'^/2kBT), where T denotes temperature and ks the Boltzmann 
constant and Ji^2 are the interaction constants. Although the model (13) looks like asym- 
metric with respect to the second subscript m, actually this asymmetry does not influence 
on the final result. It is obvious, that at Ji = the Hamiltonian (13) describes ID Ising 
model at small magnetic field. However, the situation is quite different, if J2 = and we 
have two-dimensional Ising model with the interaction constants Ji and J| = H'^/ksT. For 
this model as it is stressed above, there is no phase transition at finite temperature because 
this modell actually describes ID Ising model in small magnetic field H. It is obvious that 
the including an additional nearest neighbours interaction J2 7^ corresponding to 2D Ising 
model in a small magnetic field {H ~ e -C 1), does not lead to the phase transition also. 
The last one can be derived directly from the Hamiltonian (13). The same reasoning are 
also valid for the 3D Ising model in a small magnetic field; here it is also possible to in- 
troduce Hamiltonian, just like to (13) with the field-square term, using for this case the 
representation: 



Zs{h) = (2cosh/.)^^^(0|T3T2TiT,|0), 



Ti = exp 
T2 = exp 
T3 = exp 

and (/i = tanh^(/i)) 

r rNMK N-n 



N,M,K 
n,m,k=l 



^nmk 0inmk){0i 



t 

n+l,mk 



^NM,K ,0] 



-^2 Sn,,m,fc 

^N,M,K 



+ f3n 



nmk 



"nmk 



E E 



nmk n+s,mk 



NMK M-m 



~^ ^ y ^ y ^nmk^n' ,m+t,k ~^ ^ ^ ^ ^ ^nmk^n'm' 



nm,k+l , 



NMK K-k 



k+l 



nmk s=l nn'mk t=l nn'mm'k 1=1 

for the 3D Ising model in an external magnetic field |l^, where anmk, ■■■ are Fermi operators. 



III. CONCLUSIONS 



Despite the three representations (l)-(3), (4)-(5) and (6)-(8) for the partition function 
of the Ising model in an external magnetic field include "unpleasant" phase factors, the 



5 



representation (4)-(5) has some advantages comparing to the (l)-(3) and might be to (6)-(8) 
ones. To the author's mind, the apphcation of well-known approximate methods (analytical 
as well as numerical) to the Hamiltonian (13) can lead to the much better results than it 
took place in case of standard Hamiltonian of Ising model in an external magnetic field (with 
linear-field term). It is also possible that the application of the direct Onsager method |1|,|^ 
to the operator V = V1V2V/1 (4)-(5) could allow to find out the eigenf unctions and eigenvalues 
of that operator. 

I am grateful to I. Tralle for his assistance in preparation of the final form of this paper. 
This paper was supported by the KBN grant Nr 2 P03B 088 15. 
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